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Abstract 

Q_i' We argue from two complementary viewpoints of Holography that the 2-point correla- 

tion functions of 1/2-BPS multi-trace operators in the large- N (planar) limit are nothing 
but the (Wick-rotated) S-matrix elements of c = 1 matrix model. On the bulk side, 
we consider an Euclideanized version of the so-called bubbling geometries and show that 
the corresponding droplets reach the conformal boundary. Then the scattering matrix of 
fluctuations of the droplets gives directly the two-point correlators through the GKPW 
prescription. On the Yang-Mills side, we show that the two-point correlators of holomor- 
phic and anti-holomorphic operators are essentially equivalent with the transformation 
functions between asymptotic in- and out-states of c = 1 matrix model. Extension to 
non-planar case is also discussed. 
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1. Introduction 



A certain class of correlation functions of M = 4 super Yang-Mills theory can be com- 
puted exactly using a free-field approximation. A very interesting nontrivial example 
of such cases is the two-point correlator of multi-trace local operators belonging to the 
1/2-BPS sector. In this case, the free-field approximation of the gauge theory can fur- 
ther be mapped to a one-dimensional matrix model [U [2] . Hence, such correlators can 
be conveniently classified by using the language of one-dimensional non-relativistic free 
fermions. It is quite remarkable that the similar fermion-liquid ('droplet') picture also 
arises in the holographically dual description, namely, on the supergravity side [5] as the 
classification of a class of the so-called bubbling geometries with the same amount of (su- 
per)symmetries as the 1/2-BPS sector of the gauge theory. The correspondence of both 
sides in this particular sector can be regarded as yet another evidence for the AdS/CFT 
correspondence. 

It would be of some interest to see whether this correspondence of classification at the 
level of spectrum could be further extended to agreements of physical observables. The 
purpose of the present note is to provide such an example. According to the standard 
prescription of the AdS/CFT correspondence, the gauge-theory correlators are interpreted 
as the amplitudes of propagation of supergravity modes in the bulk, connecting conformal 
boundary to conformal boundary in the AdS background [4]. However, it has not been 
clear what should be the corresponding interpretation of the above 1/2-BPS two-point 
correlators for general multi-trace operators. 

We will argue that a reasonable holographic interpretation of two-point correlators 
naturally emerges if we consider the bubbling geometries in an Euclideanized picture. 
The two-dimensional droplets characterizing the bubbling geometries are then infinite do- 
mains which extend to the conformal boundary of the asymptotically EAdS background!)] 
As opposed to a circular droplet which is located close to the center of geometries in 
the Lorentzian case, the ground-state droplet is an infinite wedge region bounded by a 
hyperbola. This can be related to the fermi sea of the ground state of the c = 1 matrix 
model. Since the excited states of the droplet propagate along the hyperbola, we can 
naturally identify the scattering amplitudes of excitation modes along the hyperbola to 

* This possibility was first emphasized in [5] . 
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be the correlators on the gauge theory side through the standard interpretation. 

On the other hand, from the viewpoint of the c = 1 matrix model, we can argue 
that the correlators of 1/2-BPS multi-trace operators which take a form of correlation 
between holomorphic and anti-holomorphic matrix operators is reinterpreted as defining 
an S-matrix as the amplitudes of transformation between incoming and outgoing multi- 
particle states of excitations near the fermi sea. 

In the next section, we first describe an Euclideanized version of bubbling geometries. 
Then in section 3, we discuss the scattering amplitudes (which we call 'boundary S- 
matrix') of excited modes along the ground-state droplet in classical approximation. In 
section 4, we discuss some examples of matching between boundary S-matrix and the 
correlators. We also study the limit of large R-charge momentum. The latter will explain 
the origin of a previous partial proposal due to ref. [6 J for a possible relation of the 1/2-BPS 
correlators with the vertex of collective field theory of the c = 1 matrix model. In section 
5, we discuss the S-matrix interpretation of 1/2-BPS correlators from the standpoint of 
the c = 1 matrix model per se. This allows us to provide an intrinsic interpretation of the 
correlators as an S-matrix within a logic of matrix model. We also provide a nontrivial 
example of higher-genus effect which shows the correspondence of 1/2-BPS correlators 
with the c = 1 S-matrix to arbitrary genera in the large momentum limit. 

2. Euclideanized bubbling geometries 

It is useful to first recall how the ground state, AdSsxS 5 geometry, is embedded in the 
LLM metric [3] in the Lorentzian case, 



All of the functions in this metric (and also the RR-fields which are suppressed here) are 
determined § by specifying the value, either 1/2 or —1/2, of a scalar function z(xi,y) on 
the plane D at y = 0. Under this boundary condition, the Laplace equation for z(x{,y) 



ds 2 = -h~ 2 (dr + Vidx 1 ) 2 + h 2 (dy 2 + cfxW) + ye G dtt 2 3 + ye 




(2.1) 



didiz + ydyi-*-) = 



(2.2) 



y 



is solved as 




(2.3) 



§ The equations are h 



2ycoshG, z 



itanhG, yd y Vi = eijdjz, y(diVj - djVi) — Cijd y z. 
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The ground state corresponds to a circular droplet of radius r$ = y/4jrgJV 



■ -1/2 (xi.Xo) G circular disk of radius r , n . s 

*(^°)H +1/2 otherwise ^ 



which gives 



1 /" r'dr 

2(^1, »2,3/) ~ = 



2 K ./Disk r2 + r ' 2 — 2rr' cos 0' + y 2 } 2 

2 




r 2 — rl + y 2 1 



(2-5) 



2^{r 2 + r 2 + y 2 ) 2 - Ar 2 r 2 2' 

using spherical coordinates for the 2- dimensional plane (^1,^2) = r(cos0, sin0). Then 
defining new coordinates p, 9 and by 

y = r sinh p sin r = r cosh p cos 0, = — r, (2.6) 

the LLM metric above reduces to the standard AdSs x S 5 metric expressed in terms of the 
global coordinate, 



ds 2 = r 



cosh 2 pdr 2 + dp 2 + sinh 2 pdn 2 3 + dO 2 + cos 2 #d0 2 + sin 2 6dtt 2 3 . (2.7) 



Let us now perform a (double) Wick rotation [7] r — > — ir, — ► — ("0 = ^ — r — > 
— i-0) under which both the AdS metric and the RR- fields are transformed 'covariantly' 
into the Euclideanized AdS (EAdSsX S 4 ' 1 ) background with the metric, 



ds 2 = r 



cosh 2 pdr 2 + dp 2 + sinh 2 pdQ 2 3 + d6 2 - cos 2 6dip 2 + sin 2 6dQ 2 3 



X8) 



Since the signature of this metric is still 9+1 in 10-dimensional sense, supersymmetries can 
be preserved by a suitable renaming of spinor variables. The two-dimensional coordinates 
(xi,X2) are then transformed as 

X\ — > X\ = rcosh-0, X2 — > 1x2 = irsmhijj. (2-9) 

This exercise implies that for discussing generic Euclideanized LLM ansatz, it is sufficient 
to make the double Wick rotations X2 — > 1x2 and r — > — ir. The vector field Vi must also 
be rotated covariantly as V\ — > —iV\, V 2 — > V%. 

Therefore the Laplace equation now becomes a hyperbolic wave equation 

[d\ - dl)z + vdyfiZ) = 0. (2.10) 



Although it is actually sufficient to consider a wedge region xf — x% > 0, x± > which we 
hereafter denote by the symbol W, for discussing the asymptotically EAdS backgrounds, 
let us consider this wave equation in the whole (Wick-rotated) (xi,^) plane which we 
denote by the same symbol D as in the Lorentzian bubbling geometries. The region 
outside the W will be denoted by W = D— W. We have to specify the boundary condition 
at y = for the whole plane D just as in the Lorentzian case by demanding smoothness 
of solutions at y — 0. 

The solution must satisfy the boundary condition z(xi, y)\ y =o = z(xi, 0) = 1/2 or —1/2 
at y = 0. Under this boundary condition, the appropriate solution which includes the 
above EAdS background as the ground state is 



z(x 1 ,x 2 ,y) = Im 

7T 



dec -| dec o r/ 

D [{Xl 



z(x[,x' 2 ,0) 



X 



l\2 



[X2 



X 



r\2 



This is verified by noting that 

,2 



lim 



V 



+ y 2 



ir 



(2.11) 



x[) 2 - (x 2 



y^o L[(xi — x' x ) 2 — (x 2 — x' 2 ) 2 + y 2 — ie} 2 [(xi 

= —2iri5(xi — x' 1 )5(x2 — x' 2 ). 
We here chose the boundary condition outside the wedge region W as 



' ^2 



+ y 2 



iel 2 J 



(2.12) 



z(x[,x' 2} 0) = 1/2 for (xi, x' 2 ) G W. 



(2.13) 



After obtaining solutions, we can restrict ourselves to the wedge region (xi,^) £ W for 
general nonzero y, which is related to the Euclidean plane D of Lorentzian solutions by 

For example, the EAdS solution discussed above is obtained by replacing the circular 
disk (with value —1/2) of the Lorentzian case by an infinite domain H (0 < r < ro) in 
the wedge region W bounded by a hyperbola at r = r . It can be checked by explicitly 
evaluating the above integral that the expression of z for the EAdS takes the same form 
as in the Lorentzian case 

' * z(x[,x' 2 ,0) 



— Im 

7T 



dx-^dx<2 r/ 

D [[Xi 



X 



1/2 + — Im 

7T 



doc -| doc n r / 

[(*i 



' ' 2 - (x 2 - x' 2 ) 2 + y 2 
1 



x' x ) 2 — (x2 — x' 2 ) 2 + y 2 — ie] 2 



r 2 — Tq + y 2 



2 a/ (r 2 + Tq + y 2 ) 2 — Ar 2 rl 



1/2 if r > r , y = 
-1/2 if r < r , y = 



(2.14) 



as it should be, since the ground state droplet in either metric has no explicit dependence 
on the angle <fi <-» iip. 

Thus we can define Euclideanized bubbling geometries by setting hyperbolical droplet 
at y = in a similar way as in the Lorentzian case. The excited droplets can have 
arbitrary holes and islands of droplets in the wedge region W. The hyperbolical shape 
of the ground-state droplet implies that in comparing with the free-fermion picture of 
matrix models, we should consider the usual c = 1 model with inverted harmonic oscillator 
potential, instead of the ordinary harmonic potential of the Lorentzian case. There are 
two regions (corresponding to positive or negative values of X\) separated by the inverted 
potential. But for discussing deformations of the ground-state hyperbola in the classical 
approximation, it is sufficient to consider only one of them. Since after all we are discussing 
the classical supergravity solutions which are only justified in the large N limit, we are 
allowed to treat the fermi sea in the semi-classical approximation on the matrix-model 
side too. 




Fig. 1: The shaded region (H) in the Fig. 2: In the fermion picture, the ground 

wedge region W in the plane D represents state droplet corresponds to filling levels 

the ground-state droplet. from the top of the potential. 

There is actually a subtle difference that the droplet should be considered only in the 
wedge region W. In terms of fermionic picture, this amounts to considering one-particle 
levels only up to the top of the inverted harmonic potential. See Figs.l and 2. The levels 
above the potential must be totally ignored, and the ground state is the one in which 
the levels are filled from the top of potential, not from the (infinite) bottom of a suitably 
regularized potential as usual. In other words, particles and holes are interchanged there, 
comparing with case of the usual c = 1 matrix model. Clearly, however, these differences 
do not matter again in the semi-classical treatment of droplets. 

Now going back to the EAdS geometry obtained by our Wick rotation, the wedge 
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region W of the droplet plane reaches the conformal boundary corresponding to the limit 
ij} — ► ±00 which implies p — > 00. In discussing the boundary-to-boundary propagation of 
supergravity modes in the AdS background, it is usually most convenient to use Poincare 
coordinates. Therefore let us try to define the corresponding coordinates for general 
(Euclidean) LLM geometries. Since we will be interested in deformations of the ground- 
state droplet in the fixed plane W with two SO4 symmetries kept intact, it is sufficient 
to deal with four coordinates (r,xx,X2,y). They can be reparametrized as (r,p,ip,9) by 
defining 

x± = r coship, X2 = r sinh^, r = r cosh p cos 6* , y = r sinh p sin 9 (2-15) 

which are the same coordinates as we have used for the EAdS background. Then, for 
the two-dimensional part described by (r, p), we can make a coordinate transformation to 
(it, V4) defined below, where the direction of V4 is interpreted as the 4-th direction of the 
conformal boundary, by assuming that we restrict ourselves to the bubbling geometries 
which reduce to EAdS background asymptotically as u — > 00: 



In the case of the EAdS, this two-dimensional part reduces to the two-dimensional section 
(it, Ui) of the Poincare patch (u, i>i, v 2 , v 3 , U4), whose metric is nothing but the EAdS 2 , 
u~ 2 du 2 + u 2 dv\. as is easily seen by comparing this definition with the ordinary Poincare 
coordinates of the EAdSs background. 

Eqs. (12.1 5p show that in terms of the coordinates p and 9, the droplet plane D at 
y = is described by two patches, either p = or 9 = as in the case of Lorentzian 
convention. Inside the region W, the former corresponds to < x\ — x\ < Tq, while the 
later to x\ — x\ > Tq. The EAdS ground state corresponds to the droplet where the first 
patch is completely filled. The droplets of general excited states have holes inside the first 
patch, and the second patch has nonzero occupied regions. If we take into account the 
Euclidean time coordinate r, the boundary along which the two patches are sewn is the 
product space of hyperbola x\ — x\ = Tq and a semi-circle defined by 




(2.16) 
(2.17) 



cosh p sinh r 



UV4, 




(2.18) 



— + v 4 = 1 (it = coshr, V4 = tanhr). 



(2.19) 



u 
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The latter equation shows that the two-dimensional section in consideration reaches the 
conformal boundary u — > oo at v 4 = ±1. If we wish, this can be regarded again as 
a hyperbola by defining new variables q = u, p = 11V4 in terms of which the above 
equation is expressed as q 2 — p 2 = 1. We note that the semicircle coincides with the 
'tunneling trajectory' [5] which describes two-point functions of the BMN operators 
semi-classically using the GKPW relation. The time r plays the role of affine parameter 
along the trajectory. 

3. Scattering of droplets 

3.1 Boundary S-matrix 

Consider first a form of droplets which can be represented by the profile of their boundary 
in W 

x\ = a(ip) cosh-?/', x 2 = a(V0 sinh^>, (3.1) 
a(V) =r + a(VO- (3.2) 

The function a(ip) measures the deformation of the boundary profile of the ground-state 
droplet. In terms of the coordinate ip — ip — r of the EAdS background, the above droplet 
should be reinterpreted cLS cL 7~~ dependent form 

X\ = a(ip) cosh(r + ip), x 2 = a{ip) sinh(r + ijj). (3.3) 

In what follows, we rename ip by ip for notational simplicity. The rationale for this 
interpretation of time dependent droplet is that the equation of small deformations around 
the EAdS background must be 

£ = (3 - 4) 

which is obtained by our Wick rotation from the corresponding equation of AdS back- 
ground as established in ref. [H] by studying the quantization of fluctuations around the 
AdS geometry. The ground state profile does not change under this motion, while the ex- 
cited state profiles moves along the hyperbola, as is familiar in the treatment of collective 
motions of fermi liquid in the inverted harmonic oscillator potential in the semi-classical 
approximation of the c = 1 matrix model. The equation of supergravity around the EAdS 
background is obeyed by solutions obtained from these time dependent droplets at y = 
through f[2~TTD . 
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This implies that the scattering amplitude of the excited states of Euclideanized bub- 
bling geometries from conformal boundary to conformal boundary in the bulk essentially 
coincides with that of the deformed droplets propagating from r = — oo to r = +00. This 
is also consistent with the fact that the energy of the droplet in supergravity is directly 
given |3] by the collective Hamiltonian of the droplet, together with the known observa- 
tion that the symplectic structure of supergravity action also reduces [9] to that of the 
collective fields of the droplet. 

Then, according to a general idea of the holographic correspondence, the GKPW rela- 
tion, between bulk-boundary propagators and correlators on the boundary, we naturally 
expect that the S-matrix of droplets will be directly connected to the two-point func- 
tions of multi-trace 1/2-BPS operators of the (Euclideanized) SYM, under an appropriate 
mapping between the deformations of the ground-state profile and the set of those local 
operators. 

There is a well-known method (JO] for deriving the S-matrix of droplet in this picture 
which has been established in the context of c = 1 matrix model. To make the present 
paper reasonably self-contained, we briefly review the basic idea. Renaming X{ by x\ — > 
x, x 2 — > p, the equation of motion for the profile in the phase space (x,p±(x,r)) is 

d d 

— p± = x- P±g^P± ( 3 - 5 ) 

where we interpret the momentum as a function of (x, r) and the suffix + and — indicate 



two regions p + > and p_ < 0, respectively: p± = ±a/x 2 — a 2 . The S-matrix is by 
definition given by considering the relation between two asymptotic regions r — ► ±00. 
Following [10], we set for sufficiently large x 



x = e q , p ± = ±e q Te- q e ± (q,r). (3.6) 

Since in these asymptotic regions we have x ~ °W.q-t-iI> anc i ^ ^MLqT+iP f or T _^ zpoo 
respectively, the e fields behave, with respect to the dependence on r, as e±(q, r) ~ 
e±(r =F <?)• Hence the lapse of time between incoming (7$) and outgoing (7*) waves at a 
fixed large value of x = e q satisfies 

a(ib) 2 , , a(ih) 2 

Tf _ Ti = 2q + \og^L, e .(q,T t ) = e + (q,r f ) = ^L (3.7) 

which lead to a functional equation for asymptotic forms 

e+ (r -q) = e_{T-q- log ^=^). (3.8) 
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Or, using the fluctuating fields 5± defined by = c o + $±( x )i with Co being the scale 
of the static ground-state droplet, 

5+(x) = 5.{x - log(c' + 5+(x)). (3.9) 

The solution is given as (see [11] for more details) 

This can be interpreted as the relation between in- and out-fields of collective fields 
X±, defined by 

S ± ~ V^(d T =F d q ) X ±, X±(r TO) — —±= I a ± (Oe*^^- (3.11) 

V47T y 4 

satisfying the canonical commutation relations, 

[a±(£W(0] = -£<*(£' + £), o±(-o;)|0) = o;>0. (3.12) 
In terms of the normal-mode operators in the momentum representation, we have 

°*M - g (Ir'r^W / ' {,( ' -E«(nM«) (3-13) 

and the S-matrix in the classical (tree) approximation is given by 

*(E W « ~> E<4) = <°in«-(-^]>+^)|0>. (3.14) 

« i 

For example, for n — > 1 scattering (n > 2), the S-matrix elements are, up to the 
delta-function of energy conservation (u — u)\ + ■ ■ ■ ou n ) which we will always suppress in 
what follows, 

/ 2 \ n-l 

(0\a-(—u)oc + (ui) ■ ■ ■ at + (a>„)|0) =>- ( — J (—iu)(—iu - 1) ■ ■ • (—iu - n + 2) u)\ ■ ■ -u n . 

Vc / 

(3.15) 

In applying these results to our case, we have to Wick-rotate the momentum (=±energy 
on the mass-shell) as £ = u — > iJ with J being the R-charge, in order to take into 
account our prescription of Euclideanized bubbling geometries. Note that this procedure 
effectively makes the coupling constant, 1/cq, pure imaginary, apart from an overall factor 
i which can be absorbed into Wick rotation of the 5-function of energy conservation. 
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3.2 Hamiltonian formalism 

We finally recall that the same S-matrix elements as above are obtained by the Hamil- 
tonian formalism of collective field theory [12]. The equation (13.51) can be recast in the 
Hamiltonian form 

d rP± = i[H,p ± ] (3.16) 
by defining the effective Hamiltonian and the commutation relations as 



H 



dx 



[x 



pi 



(x 2 + jj)p. 



6 2 / V 6 

\p±(x,r),p±(x',r)] = TiS'(x 



x 



(3.17) 
(3.18) 



Here, /i is an arbitrary constant, corresponding to an integration constant for the equation 
of motion. Note also that we use the usual Lorentzian convention in writing down the 
equations of motion, by interpreting H3.4[) as being due to the inverted harmonic potential 
V(x) = —x 2 /2. Actually, there is an ambiguity whether we take positive /i or negative one. 
Here we choose the positive convention for definiteness. This choice has an advantage in 
that the interaction Hamiltonian is not singular. The negative fi, which looks more natural 
in view of the form of the profile function, would give a singular interaction Hamiltonian. 
Remarkably, however, there exists a duality that both give the same S-matrix with suitable 
regularization for the negative choice. 

To make the system look more like a usual canonical system, we define the shifted 
fields <f>± (the sign of /i becomes relevant here): 



(x, t) = ±(^x 2 +Jl + cf)±(x, r)), 



H 



P±{X,T 
1 



dx 



2 v / x 2 + M0+ + 2 ) + () 



1,7 



3 + 3 -) 



which reduces, by further making a change of variables x = /i sinh a, 



H 



da 



i 

+ 6 



da 



dx 



1 + <t 



(3.19) 
(3.20) 

i±, to 

(3.21) 



Here we have subtracted a (field- independent) c-number contribution. Then, using the 
normal-mode expansion in the interaction representation, 

, , 1 



/27T 



-i£(Y=F<x) 



"(0 



(3.22) 
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with [«(£), «(£')] = —uj5(!; + £') which can be identified with the normal-mode operators 
of x± introduced above in the asymptotic region \a\ ~ q — > oo, we obtain 

H = H 2 + H 3 (t), (3.23) 

/■oo 

# 2 = / d£a(Oe*(-0, (3-24) 

1 /"°° 

F 3 (r) = - y d 3 e /(6 + 6 + 6) e"^+^M6 W6)a(6) (3.25) 

with 

/(0 = r da -^^-e^ = -^f-. (3.26) 
J-oo pi cosh a /i smh 7r£ 

Note that the range of the a-integral becomes the whole real axis by combining <f>± con- 
tributions into a single integral, even though the original range was the half real axis 
< a < oo. If we choose negative /i, the form factor in H3 would have been 1/ sinh 2 a 
instead of 1/ cosh 2 a. Using this interaction Hamiltonian, it is straightforward to compute 
the S-matrix in perturbative expansion in 1/fi oc 2/c 2 ,. The agreement of the results with 
the polynomial form exhibited in (13.131) and (13.151) has been confirmed in [T3] [14] [15] . 
together with some extensions to quantum corrections. This is somewhat miraculous in 
view of the presence of a nontrivial form factor /(£) in the interaction Hamiltonian (13.251) . 

4. Comparison with the two-point functions of multi-trace operators 

Let us recall the structure of a representative set of the chiral primary 1/2-BPS opera- 
tors on the Yang-Mills side, which are characterized SO (4) symmetry and the conformal 
dimensions A = J, 

i 

where Z = (05 + icj)^) / y/2 is the complex scalar field with a unit R-charge J = 1 with re- 
spect to the rotation in the 5-6 plane. Due to the non-renormalization property, two-point 
correlation functions of these operators with their conjugate set of operators constructed 
in terms of Z = (0 5 — i(j) e )/y/2 are given by the free-field results, 

(^,J5...^(^) O&x^,...,^)^)) = ^({(^)» C- 7 )}. ^0^4(a;, a;')-' (4.2) 

where D^x, x') oc \x — x'\~ 2J with J = J { = £\ J'i is the massless free-field propa- 
gator in 4 dimensions and the function F({(J'), (J)},N) is determined by the free-field 
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contraction among the indices of scalar fields 4>i between J (x) and O (x 1 ). Obviously, 
the function F is completely independent of spacetime coordinates. Furthermore, it is 
independent of the spacetime dimensions and signature. The double Wick-rotation we 
have discussed in the previous section requires us to Wick-rotate the angle coordinate in 
the 5-6 plane. In terms of the above scalar fields, it amounts to rotating 6 to pure imag- 
inary axis, i(pQ, and hence replace Z (Z) by Z = (0 5 — 6 ) / \[2 [Z = (0 5 + 6 ) / \^2) with 
6 now being quantized with negative metric. The crucial property (Z Z) = (Z Z) = 
is preserved under this procedure, and hence this rotation does not affect the above final 
form of the correlators. 

For the specific purpose of the present section, comparison of the S-matrix elements 
with the correlators, it is sufficient to study only the function F({( J'), ( J)}, N) which is 
independent of possible choices of matrix models. Let us therefore suppress the spacetime 
coordinates and, hence, spacetime dependent factor such as D A (x,x') in what follows. 
The simplest case m — n — 1, we have 

(TrZ J TrZ J ) = JN J . (4.3) 

The results of the previous section suggests that this trivial two-point function should be 
interpreted as the trivial (1 — > 1) S-matrix element 

(0\ct-(—iJ)a+(iJ)\0) = J. (4.4) 

As before we have suppressed an overall 5-function factor —i5(J — J) = —iS(0) of mo- 
mentum conservation, which is common for all S-matrix elements. Then, the natural 
normalization for making correspondence between the S-matrix elements of the droplet 
and the correlators is a(iJ) <-> jpj2TrZ J for incoming states and a(—iJ) <-> j^TiZ 3 
for outgoing states. 

Once single-trace operators are related to single-particle states, the multi-trace oper- 
ators must be interpreted as multi-particle states. Under this interpretation, let us study 
some examples of (n — > 1) correlators in the leading planar approximation in the large N 
limit. We present the results of graphical computations for n = 2, 3, 4 (J = Y17=i ^»)- 

F(J, { J 1: J 2 }, N) planar = J.hJ 2 N~\ (4.5) 
F( J, { Ji, J 2 , J 3 }, N) planar = J^J^Wi - 1) + J1J2W2 ~ 1) 



13 



+ hhh{h - 1) + 2Ji J 2 J 3 )N- 2 = J t J 2 J 3 J{J - 1)N- 2 (4.6) 




where for n = 3 the first 3 contributions come from 'chain'-type diagrams in which the 
3 traces of the initial state connected like a chain of 3 beads, while the last one comes 
from a 'clover '-type diagram in which all of the 3 traces are connected simultaneously at 
a single point. 

oooo 9P 

Fig. 3: Four types of planar contractions of multiple traces into a single trace for 1 <-» 4 
scattering. 

F(J, {Ji, J2, J3, J4}, N) p i anar = J 2^Ji J 2 (J 2 — 1)J 3 (J 3 — 1) J4 

+ J X J3(J3 " 1)Wa - 1) J 2 + hW* " 1)Wa ~ 1)J 3 
+ JsJ^Jx - 1) J 4 (J 4 - 1) J 3 + J 3 Jl(Jl " 1) W* ~ 1)^4 

+ 6(jxJ 2 J 3 (J 3 - 1) J 4 + -h-hWl ~ 1)^4 + ^1^(J 2 - 1)^3^4 
+ J1J2J3MJ4 ~ 1)) 

+ (jiCJi - l)(Jx - 2) J 2 J 3 J 4 + ^1^2(^2 - 1)(J 2 - 2) J 3 J 4 
+ ^72^(^3 " 1)0/3 - 2) J 4 + -h-h-h-h{-h - 1)(J 4 - 2)) 
+ QJ^JsJ^N- 3 = J X J 2 J Z J±J(J - 1)(J - 2)A^~ 3 (4.7) 

where, similarly with the case n — 3, the first round bracket comes from chain diagrams 
of 4 beads, the second from three-leaf clover with one of the leafs being a chain of two 
beads, the third from 'sunflower '-type diagrams in which three traces are connected to 
one central trace at three separate points, and the last from 4-leaf clover diagrams in 
which all of 4 traces are connected at one point, respectively. See Fig. 3 for illustration. 
These results precisely match to the corresponding S-matrix elements (13.151) . provided the 
expansion parameter is related by ^2/cq —>■ l/N. 

The same results as above can also be obtained using the exact general formulas [17] , 
which are expressed as linear combinations of the ratios of Gamma functions, such as (say 
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71=1,2) 



^^ww^u _ 1 ^(iY + j 1 + j 2 + i) , r(iv + i; 



(Tr(Z )Tr(Z Jl )Tr(Z j2 )) 



J+lV T(N) r(iV-J x -J 2 ) 

r(jv + j 1 + i) r(iv + j 2 + i) x 
r(jv-j 2 ) r(iv-JO /' 1 j 

for these type of correlators. See the Appendix. These exact formulas will also be used 
for a study of the higher-genus effect in the next section. 

The computation of higher-point cases and more general configurations of initial and 
final states becomes increasingly cumbersome. To convince ourselves the validity of agree- 
ment further, it is useful to consider the limit of large momentum J. Let us define a unitary 
scattering operator S = exp V by 

n times n times 

Sa+^S- 1 = £-L'[V, [V, [■ • • , [V, a+Mj^ffl = a_fo). (4.10) 

n=0 

It should be kept in mind that here the commutators should be understood as Poisson 
brackets, since we are in the tree approximation where the operator ordering is irrelevant. 
In the limit of large momentum, we can replace the (13 . 13[) by (k = 2/c 2 ,) 

a ± (r/) = fy^i^p ! ^r/-5>)(nM&))- (4-H) 
P =i p - J 

Then, we find that the scattering operator can be written in a simple closed form as 

S = exp [i~ j d£i J d& J d& 5(£i + 6 + 6)a + (6)a+(6)a+(&)] (4.12) 
For example, the second term n = 1 (p = n + 1) is 

[*g / / ^6«+(6)a+(6)a+(-6 -6),a+(^)] = -y / ®+(£i)u+(v - CO = Oi 
The third term (n—2) is equal to 

2^6 / ^ / ^ 2a +^ 1 ) a +^ 2 ) a+ ( - ^ 1 - ^ 2 )' _ 2 _ / ^ ia +^ 1 )"+( r /-6)] 



in 2 rj 



dti I d& (^i)a+(6)«+(6)a+(^ - 6 - 6 



2j 



J d& y d6 y ^3 5(6 + 6 + 6 - r/) i 2(6 + 6 + &)«+(&)«+(&)«+(&) 
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= J <%! J d& j d&5fa + £ 2 + £3 - V) a+(6)« + (6)« + (e 3 ) = 2 

To establish the general terms in the claimed form, we can proceed by mathematical 
induction. Suppose the n-th term in the expansion is given by 

(ir)) n K n f f 



(n + 1)! 

Then by similar manipulation as for n = 2, we can check that the next term n + 1 is 
indeed given by 

(n+ 1) ^6 / ^ y ^ 2a +^ 1 ) a +^ 2 ) a +(~^ 1 - ^2),O n ] = O n+ i. 

When we use the Euclidean convention from the beginning, we obtain the same result 
by making a replacement V — > V3. 

V 3 = ^ ^(«( Ji + J 2 )a(- Ji)a(- J 2 ) + a( Ji)a( J 2 )a(- Ji - J 2 )) (4.13) 

by restricting the value of (purely imaginary) momentum by discrete integers Jj's. The 
vertex is essentially the time- integrated interaction Hamiltonian (13.251) . restricted on the 
mass-shell, of collective field theory 



V 3 ~ / pH 3 (r). 

J 27T 



In fact, it is not surprising that in the limit of large momentum, the intermediate states 
behave almost like on-mass-shell states. 

Interestingly, precisely the same formula (I4.13P has previously been proposed in [6] 
which seems to be motivated purely on a combinatorial basis that this overlap-type 3-point 
vertex is simulating the structure of splitting and joining of matrix traces]^ The suggestion 
made in the last reference was that this formula should be valid for arbitrary genus in 
the BMN limit (J,N — > 00, g 2 — J 2 /N =fixed). Our results can thus be regarded as a 
generalization of this proposal for arbitrary finite J at least in the planar approximation, 
by making clear the basis for the correspondence between the c = 1 matrix model and the 
bubbling geometry. In particular our S-matrix interpretation implies that for finite J the 
above vertex would not be appropriate for evaluation of higher then 3-point correlators. 
For that one should uses the nontrivial c=l string vertex as in |13j . 



• One of the present author (T. Y.) thanks K. Okuyama for a discussion on his work. 
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5. Droplet scattering from the viewpoint of the c = 1 matrix model 

5.1 c = 1 scattering interpretation of the correlators 

Since the function / does not depend on spacetime dimensions, we can consider a matrix 
quantum mechanics to represent the same correlators. Traditionary, we use a complex 
matrix for this purpose with special constraint called the 'lowest Landau- level condition', 
as originally discussed in pQ and followed by most literature related to this subject. 

From the reduction, it is clear that one can start with the ordinary Hermitian matrix 



model (of Lorentzian signature), J dr ^Tr 
in the following normal-order prescription, 



( dM \ 2 
I dr > 



M 2 



by considering the correlators 



k-OU v Jf) U ' :O j {Ji _jJt) m :) = F({(J'), (J)}, JVpifo O 

where -Dx(r, r') oc exp(i|r — r'\) and 



A J 



(5.1) 



Oi, : ,,(r,),, = Tr(M Jl )Tr(M j2 ) ■ ■ -Tr(M J "). 



(5.2) 



The normal product symbol : ■ • • : indicates that no contraction is allowed inside. This 
interpretation was emphasized in ref. [5] to be useful for extending discussions to more 
general 1/2-BPS operators in which the other SO (6) states than (14.11) appear on an equal 
footing. 

In our case, it is more natural to consider the model with inverted harmonic poten- 
tial of negative sign, in order to utilize the standard c = 1 matrix model to which the 
Euclideanized bubbling geometries fit well, 



S c =\ 



cir -Tr[( 



dM. 



A dr 



2 + M 2 



and consider the operators 



n± = (M ± M)/V2. 



The correlators are then given as 



F({(JO,(J)},JV)£>i(r,T 



with 



[n) = Tr^^Tr^nf ) • • -Tr^n^. 



(5.3) 



(5.4) 



(5.5) 



(5.6) 
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Note that we now have Di(t,t') oc exp(— \r — r'\). The normal-order prescription is 
automatically taken into account since (n + n + ) = (n_n_) = 0. In the language of the 
complex matrix model of ref. pQ, the matrix operators IT-|- are related to the following 
canonical decomposition 

Z=-L(A* + B), Z=-L(A + B^) 

with A, A^ being replaced by n_|_ . The lowest Landau level condition amounts to elimi- 
nating the additional canonical pair (B,B^). 

For evaluating the correlators, we can use the collective-field representation. Re- 
member that in the planar approximation of matrix models, the collective field theory 
is essentially a phase-space representation of free-fermion liquid. For our purpose, it is 
most convenient to use the coherent-state representation of the phase space in which 
n_(n + ) ~ ^(A) are regarded as generalized coordinate (z) and momentum (a), respec- 
tively. 

Tr (IF!) -> J ^- j daz J = j ^z J a(z) = a-j, 
m /ttJ\ f dz f a j J f dz a{z) J+1 

It is relevant to note that there is also a (dual) coherent state representation in which 

Trn^=/5 J = f ^z- j P(z), 

In the (analytically continued) scattering picture, the operators aj or f3j will be seen to 

coincide with in (out) fields respectively. The existence of the dual representation can be 

related to the freedom of two seemingly different choices for the sign of fi, appearing in 

the construction of the collective Hamiltonian as has already been alluded to in section 3. 

The correlator then becomes 

f dzi a(zi) J ^ +1 f dz m a(z rn ) J ™ +1 
(0| J - -jr^ - j JL + 1 ^a. j2 ■ ■ ■ a.jM. 

In this picture, the operators Tr (rr£) are simply creation operators while the operators 
Tr (n^ ) are nontrivial polynomials. It is a theorem that we shall prove below that they 
are generated by the analog of the (c = 1) ^-matrix operator. We have 

Tr (U J + ) = S^ajS. 



Consequently, the above matrix elements are 



( ) = (°I"J( ■ • • otj' m Sa-. Jl at-j 2 ■ • -a-jJO) 

To demonstrate that S is nothing but the S"-matrix of the c = 1 theory, we recall 
the analysis through collective field theory [T2]. One considers 

TrlT+ -> Tr ((P + X) /V2) ik = (5.9) 

which is represented as 



27T \z'A; + 1 z'A; + 1, 

This operator has an exact time evolution, being an eigenstate of the Hamiltonian H = 
Tr(P 2 -X 2 )/2 = Tr(n_n+ + n+n_)/2. A multiplication by e~ ikt makes it into a constant 
of motion. Evaluating this operator at early (t — > — oo) and late (t — > +oo) times gives 
the S"-matrix of the theory. One has 

a±(x) = ±x =F - (fJ. =F &±(t)) 
x 



with x ~ a/ /i/2e T . At early time, we will have only a left moving wave packet given by 
a_(i + r) and at late time we have a right moving packet defined by a + [t — r). The 
evaluation of Wik was performed in detail in [16]. At early time, with the left-moving 
packet, 



w , (v^)^ 1 fdr, kT ^ (ik + l)\ 



ik + 1 J 8tt (z'A; + 1 — \ /i 

At late time (with a right moving packet) only a single mode survives in the limit giving 



87T /i 

Consequently, one obtains a relationship between the outgoing modes and the incoming 
ones 



dr 



l+ik 



dre ikT a + (r) = - J e*" (l + ^1 (5.11) 

which is equivalent with (I3.13P obtained on the basis of a droplet picture on the bulk side. 
The c = 1 S- matrix is defined by the transformation a + = S~ 1 a_S. The identification 

&-(z) <-> at(z),ce+(z) <-> 
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then establishes the statement that the 1/2-BPS correlators coincide with the scattering 
amplitudes of the c=l theory. 

We emphasize that this holds without the so-called leg factors which appear in the 
physical interpretation of the c=l matrix model . Also the analysis given above was done 
in the tree (or planar) approximation. 

We note that this correspondence is related to the observation of [18] where an iden- 
tification between finite temperature c=l amplitudes and the normal matrix integral was 
described. In contrast to the AdS/CFT interpretation that we were concerned with in 
the present work, the correspondence discussed in [18] (see also [19]) gives a special inte- 
gration measure on the complex matrix model side such that equivalence is achieved with 
finite temperature correlation functions. 

5.2 Case of higher- genus 

Thus far we have given a basis for the correspondence between the 1/2-BPS correlators 
and the c = 1 S-matrix from both sides of bulk and boundary theories at the planar 
approximation. It is then of interest to see to what extent the correspondence will be valid 
beyond this approximation. On the bulk side, evaluating string-loop effects in the (E)AdS 
background is an unsolved problem. On the side of the gauge theory and the matrix 
models, the usual fermion representation based on the positive harmonic potential gives a 
rigorous definition of the correlators for arbitrary finite N and integer J. In principle, there 
must be a version of c = 1 fermion representation with negative harmonic potential, which 
gives the identical function F for finite N and J. It would require a special regularization 
in dealing with the negative-sign harmonic potential by appropriately taking into account 
the differences discussed in section 2. We will not elaborate on such a direction in the 
present note, since the problem is rather a matter of interpretation. 

Instead, we expect on the ground of universality that the S-matrix of the c = 1 model 
defined by the usual double scaling limit gives the right answer for large N and J. In 
the tree approximation, the large J limit was not necessary, as is reasonable since the 
momentum along the R-charge direction must be conserved and internal momenta in 
the tree approximation are fixed by external momenta. It is natural to expect that the 
discreteness of J would be washed out in the limit of large J even for internal momenta. 

In the rest of this section, we present a piece of evidence for this expectation by study- 
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ing the simplest nontrivial case, 1 — > 1 amplitude, of higher-genus effect in the leading 
large- J approximation. In the coherent state representation, one has the Hamiltonian for 
a single fermion 

h = ^a + a_ + M - (5.12) 

The wave-function can be taken as functions of a_ or (the dual picture) of a + . One has 
the simple transform from one basis to another 

if)-(a+) = —= I da_e- a+a -i) + {a_). (5.13) 



2,71 

The wave-functions both obey 

^ (a + a_ + a_a + ) ip± = (id T + /i) ip±. (5.14) 

In the coherent state representation with a_ = a a + = 8/ da , one has the equation 

d 1 \ 

a— + --njij + , k (a) = kij +jk (a). (5.15) 

From the viewpoint of the fermion phase-space (x,p) , creation-annihilation coordinates 
are the null plane coordinates a± — > x± = (p±x) / y/2. The discussion of scattering theory 
in the fermionic picture can be found in [201 EE] • 

The fermionic wave-functions ip±(x + ) obey the analogous Schrodinger eqs. The equa- 
tion is solved by 

ij + , k (x-) = A k xZ l{k -^Q(x^ + B k {-x-)^ k -^e(-x-) (5.16) 

with an analogous solution for The (Gaussian) transform then leads to the relation 
between the in and out fermion wave-function with the following reflection coefficient (see 

R{ k _ = LifCj-Cfc-M)) + 7 e -q (fc-M)l (5.i7) 

where 7 specifies the boundary conditions. 

In conformal field notation, the fermions have the expansion 

V>(z) = XV^n, (5-18) 

^ + (z) = J2^n (5-19) 
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with neZ + |, and 

The reflection coefficient is defined by 



(5.20) 



(5.21) 



where ?/' m (z) = ip-(z), ip ont (z) = ip + (z). Standard Bosonization formulas then define the 
corresponding in (Out) bosonic fields. To make contact with the previous notation we 
have 

^ ± {z)^ ± {z) = a±{z) (5.22) 

The S-matrix is then given by 



(°i IIE ( R *^)n W)l h+ n • ■ -«-Jn|0>- 



(5.23) 



k=l n 



To compare this result with the 1/2-BPS correlators, let us use the known results of 
genus expansion of two-point amplitude given in [20]. Using their notation, the contribu- 
tions up to 3 loops in the limit of large momentum q are given as 



R(q, -qh ~ -^ 5 , R(q, -9)2 



-g 9 , R(q, -q)i 
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5760" ' ^ " /d 2903040 
where we have omitted a factor (r(— |g|)) 2 fi^ multiplying each term for brevity. 
On the other hand, the exact 2-point correlator reads (see (|4.8|1 ) 



(5.24) 



G(J) 



1 



J + l 



T(N + J+1) T(N + 1) 



T(N) T(N-J)_ ' 

Using the expansion formula for the ratio of Gamma functions as discussed in the Ap- 
pendix, we find 

N J+1 A 1 f0\ \fJ\ n ( J x " n 



n=0 \ / L \ / 



Only the odd n = 2k + 1 terms are nonzero, giving the expansion 

2 ( J 



N~ J G(J) « £ 



/c=0 



JN 2fc V 2fc + 1 



J\ 2fc+1 2iV 



J 2 

— sinh — . 

J JV 



or 



A^ J G(J) « J 



1 -Z 8 1 J 12 
+ A^ 2 3!2 2 + A^ 4 5!2 4 + A^ 6 7!2 6 + 



(5.25) 
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The odd denominators give coefficients 24, 1920, 322560 agreeing with the three coeffi- 
cients in the c = 1 theory, up to an overall factor i which should be combined to the 
5- function of energy conservation as in the planar cases and also to sign factor (—1)™ for 
n loops. 

We can confirm that this agreement is valid to all orders. For this purpose, it is useful 
to recall the exact expression given in [20] for the matrix model two-point function. 

- q) = T ( - H f to f ^'**-*) - r , T ^-j ] ., ) ) (5-26) 

dfi, [ y T[^-ifi) T {-\q\ + ± - in) J J 

which gives the expansion 

R(q, -q) = (qT(-\q\)y ^ { ±- - • 1 • (| g | - l)(g 2 - |g| - 1) + • • ■ } (5.27) 

In the large q limit, one notices the identity 

d 

— R(q, -q) = fi \q\ R(q, -q) 

which follows from the fact that the factor /J 9 ' gives the dominant effect, remembering 
that we keep only the leading large q term in each order of l//i 2 expansion. Consequently 
at large q (only) we can use the formula 

R(q,-q)^^-M^-R(q,-q), (5.28) 
\q\ on 



and we have 



T(-\q\)- 2 R(q, -q) « — 



T{q + \-iii) T(\-iy) 



2|?|i I V r(|-z/i) r(-g + |- «» 



r(i + z» r(-g + | + z/i). 

Using the same expansion formula for the two terms inside the round bracket as before, 
one generates at large q the terms 



oo 

-q\ n ( g\" 



(w) 



q * — ' \n 

y n=0 



2/ V 2 



Since in the sum only the odd (n = 2k + 1) terms contribute these two contributions add 
up giving an agreement with the gauge-theory correlator G(J), under the replacement 
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— fi 2 — > N 2 with an overall factor ±i. The sign is consistent with the correspondence of 
genus-expansion parameters we found at the planar level. 

The fact that ~ N implies that the fermion levels must be filled evenly from the top 
of the potential for the ground state, quite differently from the usual double scaling limit 
in defining higher genus contributions in the traditional treatment of the c = 1 matrix 
model. This of course reflect the feature emphasized in the beginning of this subsection. 

6. Conclusion 

We have argued that the two-point correlators of 1/2-BPS multi-trace operators are in- 
terpreted as the S-matrix of the c = 1 Hermitian matrix model. The correspondence of 
both sides are fairly tight in the planar approximation. However, as we have cautioned, 
the situation is somewhat different if we consider the non-planar case, where we have 
presented evidence for the correspondence only in the limit of large momentum J — > oo. 
To establish exact correspondence for finite J and large but finite N, we need a modified 
definition of the c = 1 model, in regard to the level spacings in filling fermion states 
from the top of the potential, not from the bottom, as we have alluded to in section 2. 
It must be defined such that the 'ground state' of the model becomes equivalent with 
that of the matrix model with positive harmonic potential and consequently we can have 
a one-to-one mapping even for excited states (again from top to down) to those of the 
positive harmonic potential. 

We may also hope that, in the large J limit, we can discuss the genus corrections from 
the viewpoint of pp-wave string field theory. However, at the present state of development, 
even the question of what is the right string-field vertex for this purpose is not settled. 
For some recent works related to this question, see [23] . 

We have emphasized that unlike the known correspondence between c = 1 matrix 
model and the S-matrix of 2D string theory, the leg factor is not necessary. In ref. [8], 
it has been established that we need a special leg-like factor in order to relate the 3- 
point correlator of the BMN operator with the Euclidean S-matrix in the AdS bulk from 
boundary to boundary. However, for the extremal correlators such as 

(Tr(Z Jl+j2 )(x)Tr(Z Jl )(y)Tr(Z j2 )(y')), 
the corresponding 3-point (on-shell) vertices in the bulk vanish and simultaneously the 
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leg factor diverges, so that the product of 3-point vertices and the leg factor gives finite 
correlators. In the limit y — ► y', these extremal 3-point correlators reduce smoothly to 
two-point correlators which are dealt with in the present work. Similar arguments apply 
also to higher-point extremal correlators. It is quite significant that the droplet picture 
provides us a direct correspondence for the correlators and Euclidean S-matrix in the bulk 
without such complications. It would be very interesting if there is any extension of our 
results to non-BPS operators. 

Indeed, the c=l interpretation might be of relevance to the full AdS string where 
analogous structures and c = 1 type scaling [2U [25] has been observed recently. The 
latter correspondence regarding the world-sheet S-matrix points towards capturing the 
mixing among excited but single-body string states. In contrast to this, in our case, the 
c = 1 matrix model is relevant in understanding the multi-body states of strings, restricted 
to the ground state of spin chain states. 

We noted parallels with the relation between the so-called normal (complex) matrix 
model and the c = 1 model at finite temperature [IB] . It would be of interest to clarify 
the connection further. 

Finally, we mention also the recent investigation of [26J which discusses correlators 
and topology changes in the 1/2 BPS sector. 
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A Genus expansion from exact formulas 



Here we present some formulas which are useful for studying genus expansion on the basis 
of exact formulas for correlators. First we need an expansion formula [27] for the ratio of 
Gamma functions: 

TT^ = E l r? a ~/ + 1 l,^ "M«~'-- 

r(z + 3) ' n\r{a — 3 — n + 1) 

v ' n=0 v ' 

where the coefficient function Bn \x) is the generalized Bernoulli polynomial [28] ex- 
pressed as 

k 



*m-eac)C ,+ ;~ i )»- i >'( 



x F[fc-n,fc-//;2fc + l;— — ]. (A.2) 

x + j 

In our case, the explicit expressions of the coefficients are 

B^ +l \a) = l -{a + B-l) ) (A.3) 

Bt 0+1 \c) = \(a + 3y-^a-^3 + I, (A.4) 

sJ-^Ca) = i(a + (3f - i(a + | )(a + /?) + |a + -/?. (A.5) 
In particular for a. — [3 ^> 1, we have an asymptotic form 

r(. + /3) 2^ n J 2" ' ^ 

which can easily be understood by a 'dilute-gas' approximation for the power-series ex- 
pansion of the above ratio in 1/z. Note that in this limit (a + /3)/2 is the average value 
of the coefficients of 1/z in the product n^Zg(l + f )■ This asymptotic formula is used in 
section 5. 

Using the exact formulas p2] such as (14.81) . (14.91) and, for n = 3, 

- jir( r(iV + Jl rm J3 + ~ - r -^m^r - 

r(iV + J x + J 3 +l) T(N + J x + J 2 + l) r(JV+Ji + l) 



T(N-J 2 ) T(N-J 3 ) T(N-J 2 -J 3 ) 
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I r(iV + J 2 + l) T(N + J 3 + 1) T(N+1) \ 

+ T(N- J,- J 3 ) + T(N- J x - J 2 ) T(N — Ji — J 2 — J 3 )y ' K ' } 

together with the above expansion formula flA.ll) . we can check the results of section 4, 

although the algebra becomes increasingly tedious for larger n. 
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